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Spectral  lines  have  long  been  used  as  a  density  diagnostic  in  dense 
plasmas.     Early  theories   produced  hydrogenic  spectral  lines  that  were 
symmetric.     At   higher   plasma  densities   diagnostics  could  be  improved  if 
asymmetries  were  included.     In  this  work,  several  sources  of  line 
asymmetry  are  considered.     In  Chapter  One,  the  ion  quadrupole  interaction 
is  studied.     This  is  the  interaction  between  the  radiator   quadrupole 
tensor  and  the  electric  field  gradients  of  the  perturbing  ions.     The 
field  gradients  are  calculated  in  an  Independent  Perturber  Model.     Two 
other  asymmetry  sources,  fine  structure  splitting  and  quadratic  Stark 
effect,  are  considered  in  Chapter  Two.     Results  of  the  combination  of 
these  effect   follow  in  Chapter  Three.     Non-hydrogenic  spectral  lines  are 
also  asymmetric.     In  Chapter  Four,  a  procedure  for  calculating  the 
spectal  lines  from  Helium'like  ions   is  outlined. 
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CHAPTER   I 
lON^QUADRUPOLE   INTERACTION 

A.      Formalism  with  lon-'Quadrupole  Interaction 

In  this  chapter,  we  will   consider  the  effect  on  spectral   line  shapes 

of  the  interaction  of  the  radiator   quadrupole  tensor  with  the  electric 

field  gradient   due  to  perturbing  ions.     The  form  of  this  interaction  is 

derived  naturally  from  a  raultipole  expansion  of  the  radiator  as  a 

localized  charge  distribution  p(x)   in  an  external   potential  ^(x)   caused 

by  the  perturbing  ions    (Jackson   1972).     The  electrostatic  energy  of  the 

system  is 

W  =  /p(x)$(x)d^x.  (I.A.I) 

The  potential  can  be  expanded  in  a  Taylor  series  with  the  origin 
taken  as  the  radiator  center  of  mass. 

2 

*(x)  =  $(0)  +  x.V<I>(0)  +  1  I     x.x.  |-^  *    ...  (I. A. 2) 

2  ^^   1  J  dx.dx. 

Div(E(0)  =  ^-nlz.e&ix.) ,   where  the  sum  is  over  all  perturbing  ions  and  z^e 

i      ^ 
is  the  charge  of  ion  "i".  Since  Coulombic  repulsion  will  keep  the 

perturbing  ions  from  reaching  the  origin,   (1/6)r  Div(E(0) )  can  be 

subtracted,  giving 


?     3E. 
.  *(x)  =  $(0)  -  x.E(O)  -  -^l      (3x.x.  -  r6,.)   ^(0)  +  ...         (I. A. 3) 

ij    ^  ^  'J    J 

or, 

'•  «  ■   ■'; " 

'-■'-.''  -1-  ■       "- 
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3E. 


where 


q  is  the  total   charge  of  the  radiator, 
d  =  <ip|'er|i|j>  is  the  dipole  moment   of  the  radiator  and 
Q.  .   =  <4j|-e(3x.x.   '  r^6.  .)U>   is  the  quadrupole  tensor  of 
the  radiator. 

This  multipole  expansion  is  basically  an  expansion  in  the  atomic 
radius  divided  by  the  average  inter-ion  spacing,   or   (n  aQ/Z)/R.Q  where  n 
is  the  principle  quantum  number  of  the  radiator,  a^  is  the  Bohr  radius.   Z 
is  the  radiator  nuclear  charge  and  R^o  is  the  average  distance  between 
ions    (Sholin  1969).     The  dipole  interaction  is  second  order  in  this 
parameter,   the   quadrupole  interaction  is  third  order,  etc.     For  a  Ly-6 
line  emitted  from  an  Ar^^"^  plasma  with  an  electron  density  of   10^     cm'    . 
this  smallness   parameter  is    .15. 

To  understand  qualitatively  the  asymmetry  due  to  this  effect, 
consider  the  perturbation  on  the  n-2  levels  of  a  hydrogen  atom  located  at 
the  origin,     by  a  single  ion  of  charge  +1   located  at  -z.     If  we  initially 
consider  the  perturbation  upon  the  isolated  hydrogenic  energy  levels  to 
be  due  only  to  the  electric  field  produced  by  the  ion,  the  new  states 
are,  to  first  order  in  perturbation  theory,  the  well-^known  Stark  states 
with  a  perturbing  field  E-e/z^.     Now,   if  we  take  these  as  our  new  basis 
states  and  include  the  effect  of  the  field  gradients  due  to  the        . 


-J- 


perturbating  ion,  the  result   is  shown  in  Figure   1.     The  higher  energy 
levels,  which  give  rise  to  the  blue  wing  intensity,    (uj  >  w-,  where  huy 
is  the  energy  difference  between  the  unperturbed  initial  and  final   state) 
are  brought   in  closer  while  the  lower  energy  levels,   corresponding  to  the 
red  wing   (03  <  o)-),     are  spread  further  out.     These  shifts,  when  averaged 
over  all   possible  perturber  positions,  will  lead  to  a  blue  asymmetry 
for  Lya,   at  least  out  to  the  wings.     Similiar  arguments  can  be  made  for 
higher  Lyman  series  lines. 

Now  that   we  have  some  idea  of  what  to  expect   from  the  quadrupole 
interaction,  we  proceed  with  the  general  theoretical  formulation   (Griem 
197^*). 

The  power  spectr'um  emitted  by  one  type  of  ion  in  a  plasma  can  be 
written  in  terms  of  the  ensemble-^averaged  radiation  emitted  by  one  ion  of 
that  type  as 


P(a))  =  -!iii^  I   |<a|d|b>|^  p  6(u)«Lo   )  =  ^  Koo),  (I. A. 5) 

3c^  ab  3c^ 

where  p  is  the  probability  that  the  system  is  in  state  |a>,  d  is  the 

3 

radiator  dipole  moment,  and  oj  .    -   (E  *£.  )/h.       The  second  equation 

ab  a     b 

defines   the  line  shape  function  Koj).       Using  the  integral  definition  for 
the  Dirac  delta  function,  we  obtain. 


liui)   ^  i-  I      |<aldlb>r  p  e  ^°     dt.  (I. A. 6) 

2^  i«     ab  ^ 


Note  that  the  integrand  for  negative  values  of  t  is  equal   to  the  complex 
conjugate  of  the  integrand  for   positive  values  of  t.     This  enables  us  to 


FIGURE    1 

Energy  levels  of  the  n=2  states  of  hydrogenic  ions  in  the 
presence  o"  a  uniform  electric  field  and  in  the  presence  of  a 
field  gradient. 
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write  Ko))  as  an  integral  from  zero  to  infinity: 


m 


I((o)  =  ^  Re  /  I   |<a|d|b>|%„e''"  ^^'   dt 


0  ab 


i  (a)-(jo_j,)t 
a' 


^  Re  J   I   e^''^<b|d|a>.<a|pe~^^^d|b>dt 


0  ab 


^Re  /  dt  e^-^Vc^pe-^^^d).  (I.A.7) 

0 

L  and  p  are  the  Liouville  and  equilibrium  density  matrix  operators  for 
the  radiator-plasma  system,  respectively.  The  trace  is  over  states  of 
the  radiator  and  plasma. 

Thus  far,  Doppler  broadening  has  been  neglected.  We  will  include 
this  approximately  at  the  end  of  the  calculation  by  convoluting  the  Stark 
line  shape  with  a  Doppler  profile  derived  using  a  Maxwell  velocity 
distribution.  This  approximation  assumes  that  the  change  in  momentum  of 
the  radiator,  due  to  the  emission  of  a  photon,  is  negligible  and 
therefore  the  Doppler  and  Stark  broadening  are  independent. 

Rather  than  immediately  factoring  the  density  matrix,  we  follow  the 

general  formalism  of  Iglesias  (1984).  This  will  show,  in  a  more  natural 

way,  how  the  electron  screening  of  ion  interactions  comes  about.  We 

consider  the  full-Coulomb  radiator-plasma  interaction,  V  .  which  is 

ap' 

given  by  a  sum  of  pair-wise  additive  terms.     Furthermore,  we  will  treat  a 
two  component  plasma,  which  consists  of  electrons  and  only  one  type  of 
Ion,  hydrogenic  ions  of  nuclear  charge  Ze. 

The  Coulomb  potential  between  the  radiator  and  the  perturbing  ion 
can  then  be  written  as 


N  N   „  „  2      2 

0  Q        L     LB  L    e 


where  o  signifies  perturber  type  (in  this  case,  electrons  and  one  type  of 

ion);  N  and  Z  e  are  the  number  and  charge  of  perturbers  of  type  o 
a      0 

respectively.  Expanding  this  expression  in  Legendre  polynomials,  we 
find, 


»   V  ,1 


v^a,j)=-Z^e  I  [-^--i-6^o^P^(cos9^.).  (I.A.9) 


£=0  r>       j^ 
where 


V,   =  smaller  of  Ir  I  and  Ir.l 
<  'a'     '  J  ' 


r  =  larger  of  |r  |  and  |r.| 


e  .  »=  angle  between  r  and  r . . 


This  expression  may  be  regrouped  as 


v^Ca.j)  «  ^^''(j)  +  v^°(a,j)  +  v^^Ca.j),  (I. A. 10) 


where 


<t>^  ( j )  is  the  monopole  term  given  by  Z  (Z-1)e  /r. 


•o- 


'i-- 


i 

2  r   a 


v/(a.J)  =-z/  I  ^^P^(cose^.). 


and 


il=1  r. 

"^0 


^         0 

2  V   r  ^ 


_Z  e^  y  [— ^—  -  ^  „  Jp„(cos9  .) 


r .  ^  r 


°   £=0  r  ^^'        r.  "^'   '' 


J„    a 


v/(a.j)  =  ^  a      J 


0 


r .   >  r 


The  sum  of  v^"  +  $q°  is  identical  to  v°  in  the  region  r.    >  r^,  the 
no  penetration  region,  bat  extends  this  form  to  all  values  of  r^  .  In 
this  form,  v  °  may  be  written  as  a  "product"  of  an  operator  on  perturber 
coordinates  times  an  operator  on  radiator  coordinates,  in  the  following 


manner , 


v^°(a,j)  =  M{a)*°(j) 


K=1  ^  ^ 


1;.°(^^   +  1  T  o  a  r  °(])   +  ...).         (I. A. 11) 


C^-^'(i)^-e     I     QmnVn^J^  '"'^ 


mn 


Here  M(a)  stands  for  the  set  {uj^Ca)}  which  depends  only  on  the  radiator 
state.  For  example,  y^  -  d,  which  is  the  radiator  dipole  moment,  and 
y  -  Q/6,  which  is  the  radiator  quadrupole  moment.  The  perturber 
operator,  *°(J),  depends  only  on  the  jth  pert'urber  of  type  o.   The 
term  *°(J)  represents  the  set  l4)^°(j)},  where  ^^''(j)  -  e°(j)  is  the 
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b.-:.: 


electric  field  at  the  radiator  produced  by  the  jth  perturber  of 

type  0,  4ip°(J)  -  8  e  °(j)  is  the  electric  field  gradient  at  the  radiator 

produced  by  the  Jth  perturber  of  type  a,  and  so  forth.  By  design,  \i^ 


couples  to  ((>  . 


We  next  generalize  the  usual  definition  for  an  electric  microfield 

distribution,  W(e)  =  Tr  p(e-Z  E  ),  to  include  the  field  gradients  and 

i 
higher  order  derivatives  which  are  used  in  equation  (I.A.1 1): 


V{^)    =  Tr  p6(4'-*  ),  (I. A. 12) 


where  {('i'-*  )  =  n  &(.'¥-i'^   ). 
K-1    ^     ^ 

The  "H     and  $   are  analagous  to  4)^  in  equation  (I.A.1 1).  That  is, 

K       K  K 

they  denote  the  value  of  the  electric  field,   field  gradient  tensor,  etc. 
at  the  radiator,   produced  by  the  perturbers.     In  the  above,   i>     is  an,   as 
yet,  arbitrary  function  of  only  ion  coordinates  and  f  is  a  c-number .     Now 
equation   (I. A. 7)   can  be  written  as 


Kco)   -  -  Re    f     dt  e^       Tr{d«Tr   pjdy6(4'-*   )d(-t)} 
'^  0  a        p 


fdY  W('P)  -^e   I     dt  e^'^^  Tr{d-Tr   p6('?-$   )/W('}')3(-t) } 


a         p 


-   JdY  W(4')J(a),'}'),  (I. A. 13) 


where 
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J(ai,4')    -  -  Re    f     dt  e^*^^  Tr{d-Tr   p('i')d(-t) ) , 
0  a         p 


and 


pW    =   p6(>l'-$*)/W('J').  (I.A.IiJ) 


We  may  rewrite  JCu.^)   as 


CD 

JU,^)    -  -  Re   I     dt  e^"'^  Tr  d-f  (a,4')D(t) ,  (I. A. 15) 

^  0  a 


where 


D(t)   =  <d(-t)>, 


<X>  =  f  \a,'V)   Tr  p(4')X  for  any  X,  and 


f(a,1')  "  Tr  p(?). 

P 
The  Hamiltonian  associated  with  the  Liouville  operator  can  be 

written  as  the  sum  of  a  free  atom  Hamiltonian,  an  isolated  plasma 

Hamiltonian,  and  an  interaction  term.  Noting  equations  (I. A.8)-(I. A. 1 1) 

for  the  interaction,  we  have 

N 

0 

H  -  H  *  H^  +  I   I  v°(a,j) 
^    P   0  J  =  l 

N 
-H  *H  +1  I%.^j)  +  V  °(a,j)  *  M(a)*''(j).         (I.A.16) 
^    P   0  j-1 

The  atomic  Hamiltonian,  H-,  includes  the  center  of  mass  motion  of  the 


atom  and  the  internal  degrees  of  freedom.  The  plasma  Hamilton! an,  Hp, 
includes  the  kinetic  energy  of  the  perturbing  electrons  and  ions  as  well 
as  their  Coulombic  interactions. 

We  now  make  the  static  ion  approximation,  which  can  be  expressed  as 


e'^^^pCi-)  =  pC^-) 


This  property,  stationarity ,  implies  that  the  kinetic  energy  of  the 
perturbing  ions  plays  no  role  in  the  line  shape  problem  and  may  be 
integrated  out.  Therefore,  the  delta  function  in  p(t)  allows  us  to  make 
the  replacement  $  =  $  +  4'-*  inH,  giving 

N 


0 


» 


H  =  H  +  M(a)4'  -  \   -  II   ^.'(J)  ^  I   I   v/(a.j)  ^  M(a)(I  I  *°(j  )-$") 

^  ^   0  j   "  0  j  "^            0  j  =  1 

H  +  H  +  V-  +  V,($*),  (I. A. 17) 
a    p    2    1 


where 


and 


H  =  H  +  M(a)4', 
a    a 


■^   0  J 


V2  -  I  I  v/Ca.j), 
0  j 


V^(*  )  -  M(a)(I  I  *°(j)  -  *  ) 
0  J 


f»». 
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We  have  combined  H,  and  MCa)*  because  these   involve  only  atomic 
coordinates.     Similiarly,   the  monopole  term,   I  I  <t>o   ^J ) '   '^^^^  "°^  depend 
on  internal  atomic  coordinates  and  therefore  is  combined  with  Hp. 

Let  L(a,H'),       L   ,    L^,   and  L^($   )   be  the  Liouville  operators  corresponding 

_      _  » 

to  H   ,    H   ,    V^,      and     V, (*   )     respectively.      Further  let 
a       p        2  1 

L  ($   )    =  L   (0   )    +  L       be  the  Liouville  operator   corresponding  to  the 
radiator-plasma  interaction.     Now  use  the   Zwanzig  projection  operator 
technique  to  derive  a  resolvent  expression  for     J(a),4').       The  details  are 
in  Appendix  A  and  the  result   is 


JCaj.H-)    =  -TT"''lm  Tr{d-f (a,4')[a)-L(a.<l')-B(4')-^  U,^)]     d},  (I. A. 18) 


where 


BCy)   =  <Lj($   )>, 


^  (o),'!')    =  f    \a,^)Tr  Lj($*)  pCij;)  (w-QLQ)    ^QL^C*   ). 

p  •  ,  ' 

and 

Q  =   1-P,  where  P  is  a  projection  operator  given  by  -.       .  ^ 

P(...)    -   <(...)>. 

We  wish  to  cast  the  expression  for   Koj)   in  a  more  familiar  form,  one 
which  uses  the  standard  microfield.     To  this  end,  we  must  separate  out 
the  electric  field  dependence  from  t.       Let  "¥  «=   {e,4"},   that  is,   ¥'   is 
all  of  H"     except  for   e.       Now, 
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¥(4-)  =  Q(e)W(£|'y')  (I. A. 19) 

where  Q(e)  is  the  usual  raicrofield  (Tighe  1977;  Hooper  1968)  and 
WCeIt)  is  the  conditional  probability  density  function  for  finding 
4"  given  e.   Then, 

Ko))  =  jit   Q(e)J(a),t)  (I. A. 20) 

which  is  the  desired  form,  but  now 

J(a),e)  =  JdS"  W(el4")J(a),4').  (I. A. 21) 

k 

At  this  point,  some  simplifications  must  be  made  in  order  to  arrive 

at  a  calculable  scheme. 

1.  Neglect  ion  penetration.  With  this  approximation,  the  ion 
interaction  can  be  treated  entirely  in  the  microfield  fashion 
above.  Error  is  only  introduced  for  configurations  in  which  a 
perturbing  ion  is  closer  to  the  radiator  than  the  radiator 
electron  is.  This  is  extremely  rare  due  to  the  large  Coulombic 
repulsion  of  these  highly  charged  ions. 

2.  Choose  *  ,  which  has  been  arbitrary  up  to  this  point,  to  force 

<L  (*  )>  =  0.   This  implies  that 


0  =  Tr  p(\i/)(*  +  *^  -  <I>  ) 
P 


Tr  6(tl)-*  ){Tr  p($  +  *  )  -  p.*  )/W(i|)),  (I. A. 22) 

i        e 


where  p  =  Tr  p . 


To  satisfy  the  above,  it  is  sufficient  to  define 


*         —  1  ©        i 

*  sp.   Trp(*  +*) 
e 


$^  +  p.  ''  Tr  p$®.  (I. A. 23) 

e 


4 
Obviously,   $     is  equal   to  the  pure  ion  term  plus  a  shielding  term 

due  to  the  electrons.     This  is  not  an  approximation. 


3.     Assume  that  the  entire  effect  of  the  electron-ion  interaction  is 
to  screen  the  ion  interaction,   *^ ,   in  the  manner  above.     With 

this  approximation,   the  ion  parts  of     V^(4>   )   cancel,   leaving 

«  e 

V  (*  )    =  M(a)     I     il>^(j).       To  lowest  order  in  the  plasma 

^  j  =  1 

parameter,  this  screening  is  given  by  the  Debye-Huckel 

interaction   (Dufty  and  Iglesias   1983)  and  that  form  will  be 

assumed  here.      (Even  if  the  Debye-Huckel  form  were  not  used,   the 

procedure  would  still  be  applicable  as  long  as  the  potential  is 

«  « 

the  sum  of  single  particle  terms.)     Note  that     V^C*   )   and  L^C*   ) 

no  longer  depend  on  $   ;   therefore  they  will  henceforth  be  denoted 
simply  V^   and  Lj,  respectively.     Neglect  of  electron-ion 
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interactions  also  allows  the  density  matrix  to  be  factored. 
Thus  f(a,'i')   can  be  reduced  to 


Tr   p.6{\li-i>   )Tr   p 
1  ^ae 

f(a,,|,)    =^ ^ 


Tr   p.6{\l>-i'   )Tr    p 
,        1  ^ae 

1  ae 


Tr   p 


ae 


Tr   p 
ae 


-  f(a). 


ae 


(I.A.2i4) 


'4.     For  the  plasma  conditions  and  Lyman  lines  which  we  are 

considering,  most  of  the  line  is  within  the  electron  plasma 
frequency;    therefore  the  electron  broadening  is  primarily  due  to 
weak  electron  collisions  and  a  treatment  of  ^(oj,*)   which  is 
second  order  in  Lj   is  appropriate   (Smith,   Cooper,  and  Vidal 
1969).     This  is  accomplished  by  replacing  the  L  in  ^(oj,*)   by 
Lg.     Further  simplification  concerning  the  projection  operators 
is  given  in  Appendix  A,   with  the  result 


'U,^)    -  r~\a)Tr  Ljp(t|;)(u)-LQ)"^Lj    . 
P 


(I. A. 25) 


Note  that  the  Lq  in  the  last  equation  has  an  ion  part,  however 
this  ion  Liouville  operator  acts  only  on  functions  of  radiator 
and  electron  coordinates  and  therefore  contributes  nothing.  The 
trace  over  ions  is  now  trivial,  leaving 


if- 


I  -y 
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1    «  I 


.,    9^{ui,^\>)    '  f'^a)Tr   LjP^g(a)-L(a,iJ>)-Lg)      L^ 


f   \a)Tr   I   LAi)pii^-L{a,^^>)   -  I  L   (k))        I  L   (O. 

-^  (I. A. 26) 


where  the  Liouville  operators  have  been  written  as  sums  of  single 
electron  operators.     Terms  with     J-^IL  are  zero  by  angular  average 
and  terms  with     kn  are  zero  because  L   (k)   is  zero  unless  acting 
on  functions  of  electron  k. 

^UA)    =  f"Va)n  Tr   L   (a,  1  )f  (a,  1 )  (a)-L(a,ij;)   -  L^CD)      Lj(a,l). 

1 

where     nf(a,1)   »  Tr  p^^   .  (I. A. 27) 

^2'*'^N 

5.  The  No  Lower  State  Broadening  Ass'umption  states  that  there  is  no 
broadening  of  the  final  radiator  state  by  perturbing  electrons. 
This  is  a  very  good  approximation  for  Lyman  series  lines  because 
the  final  state  (the  ground  state)  has  no  dipole  interaction  with 
the  plasma. 

6.  The  No  Quenching  Approximation  assumes  that  electron  collions  do 
not  cause  the  radiator  to  make  non-radiative  transitions  between 
states  of  different  principal  quantum  number.  This  approximation 
is  best  for  low  level  Lyman  series  lines  since  the  levels  are 
well  separated.  Although  approximations  (5)  and  (6)  are  not 
necessary,  they  are  made  here  for  the  sake  of  numerical  ease. 
These  reduce  ^  (u)  to  (Tighe  1977) 
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e/   ^      -in    f     V   iAojt,,    -iH(1)t/fi 
A^(a))..,  -  -|^  Tr  J   dt  I  e  "  V   ,e 

^^    h   1   0    i" 


9  V,„,,ei"^'^'^^f(1).  (I.A.28) 


i"i 


where  the  trace  is  over  states  of  a  single  electron,  H(1)  is  the 
Hamiltonian  for  that  single  electron  including  its  kinetic  energy 
and  monopole  interaction  with  the  radiator,  Aid  =  u  -  (Eg-E^)/-fi 
and  subscripts  i  (f)  stand  for  initial  (final)  radiator  states. 

7.  Equation  (I. A. 27)  includes  the  ion  shift  of  radiator  states 

in  ^(uj.'i')  through  LCa.H').   Dufty  and  Boercker  (1976)  found  that 
%P  (a),"?)  is  fairly  constant  within  the  electron  plasma 
frequency.   For  this  reason,  we  will  neglect  the  ion  shift  in  the 
electron  broadening  operator  by  replacing  L(a,'l')  by  L(a). 

8.  The  real  part  of  ^(to)  is  called  the  dynamic  shift  operator.  It 
has  been  calculated  by  Woltz  (1982)  and  found  to  give  a  very 
small  red  shift.  It  will  be  ignored  here. 

Most  of  these  approximations  are  standard  in  calculable  line 
broadening  theories.  In  fact,  if  we  were  now  to  neglect  4"  completely, 
replacing  L(a,'i')  by  L(a,e)  and  Jdo.^")  by  J(a),E),  this  reduces  to  the 
full  Coulomb  (electron  interactions)  of  Woltz  et  al .  (1982).  Instead, 
limit  y  to  the  quadrupole  interaction.  Higher  order  terms  should  be 
successively  smaller  in  magnitude  (Sholin  1969)  and  the  next  higher  term, 
the  octapole  term,  is  a  symmetric  effect  (Appendix  B.)  to  first  order  in 


fe':' 
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perturbation  theory  and  therefore  should  not  have  a  large  effect  on  the 
line  shape.  This  yields 


Ko))  »  [dE  dE  de  dE  ds  dt     de  QCe) 
■'     XX  yy  zz  xy  xz  yz 


WCeIe   ,e   , .  . .)  J(a),e  ,e   ,e   ,...),  (I. A. 29) 

'  XX  yy  XX  yy 


We  now  make  the  following  approximation 


WCeIe   ,e   ,...)  ■►  6(e   -<e   >  ,  £   -<E   >  ,...).  (I. A. 30) 

'  xx'  yy         XX   xx  e   yy   yy  e 


That  is,  the  field  gradients  are  replaced  by  their  constrained 
averages.  This  form  is  exact  in  the  nearest  neighbor  limit.  It  is,  in 

general,  exact  to  linear  order  in  <e. .>  ,  but  introduces  second  order 

2  2 

errors  which  are  proportional  to  <e. .  >  -  <e.  .>  .   These  should  be 

^  IJ   E       IJ  E 

2 

three  orders  higher  in  the  smallness  parameter  n  aQ/Z/RQ^ . 

This  yields 

Ko))  -  JdE  Q(E)J(a),E), 
with 

J(a),E)  -  -iT~^Im  Tr{d«f(a)(a)-L(a, £,<£..>  )  -^(  a)))~^d}.         (I. A. 31) 
a  -^  ^ 
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The  Hamlltonian  corresponding  to  L(a,e ,<e^  .>^)  is 

H  -D«£  -  y  Q  <£  >  /6  where  p  is  the  radiator  dipole  operator  and  Q  is 

the  radiator  quadrupole  operator. 

By  taking  advantage  of  the  fact  that  the  ground  state  has  no  dipole 
or  qaadrupole  moments,  we  can  rewrite  the  expression  for  J(u),e)  as 


—  1 

J(u),e)    =  -  ^  Im  Tr  dIoj-CH^-ojq)   -  p-c   "     I  Q.j<e.j>/6   -  ^(oj))      . 

(I. A. 32) 


where  Wq  is  the  radiator  ground  state  energy  and  D  is  given  by  d-d,   and 
the  radiator  dipole  operator,   d,     is  restricted  to  have  nonzero  matrix 
elements  only  between  initial   and  final  states  of  the  Lyman  line  to  be 
calculated. 

The  Debye-Huckel   interaction  was   chosen  for  *   .       Therefore, 

5  E     .   (Ze)  H  e-'M^  -  ^  -  ^)  "•«•") 

^  ^  r  Xr         X  r       r^ 


and 


z  z  p3       ^^^  ^^'^       ;^'=r.       r^ 

where  y  -  cos (9)  and  X  is  the  screening  length.  All  off-diagonal  terras 
are  of  the  form  of  equation  (I. A. 33)  and  all  diagonal  terms  are  of  the 
form  of  equation  (I.A.3'4).  Since  only  the  z  direction  is  special 
a  is  along  2),  <e,y>,  -  <e,,>,  -  <ey,>,  -  0.   Therefore. 


^•^K 


"v^---        •      -■  .  -.-   -..r 


-zu- 


(0    +  Q   )<£   >   +  Q   <E   > 

XX    yy   XX  e    zz  zz  e 


-Q  <e  >  +  Q  <e  > 
zz  XX  E    zz  zz  c 


1  Q   (<e   >  -  ^<V-e>  ).     '  (I. A. 35) 

2  zz   zz  e   3     e 


For  the  Debye-Huckel  form, 

^1    '  r        • ;  _   ■ 

'  '  ■  ■  '  ■■  ,  '  ' 
Ve  =  4  ^^^^  e''^^Vr.  (I.A.36) 

X 


So 


1  „  .     ,....e-^''  ..   ^_^L^M_,  2. 


<e   -  ^  V.e>  -  (Ze)<^-5—  (1  ^  7  +  ^)(1-3p  )>^ 
zz   3     e        „3       A   „,^        t 


r-'  3X 


<f>   . 

c 


Now  J(a),E)  may  be  written 


(I. A. 37) 


J(a).e)  -  -  ^  Im  Tr  D{a.-(H^-a)o)  -  p-e  -  Q2Z^f>e/^  -^(oi)}'"'.     (I. A. 38) 
r 

The  quadrupole  moment  Q^^   is  easily  calculated  as  in  Appendix  C.  All 
that  is  needed  now  is  to  calculate  the  constrained  plasma  average  <f>^. 


B.  Calculation  of  the  Constrained  Plasma  Average 
From  equation  (I. A. 33),  the  constrained  plasma  average  can  be 
written 


<f>   =  <E    "   -r   V't> 

e     zz   3     e 


Obviously,  f  can  be  written  as  a  sum  of  single  ion  terms, 

f  =  Ze  I  F(r)(1-3y^),  (I.B.2) 

ion 

where 


w  ^   e"^/^  ,,   r   1  r^ 
F(r)  =  r-  (1  +  7  +  o  -p) 

r^      ^   ^  A 


Therefore,  we  have 

Jd^^r  (I  f^)6(e-E)e-^^ 


Njd3Nr  f^6(M)e-S^ 
Z^Q(e) 


(I.e. 3) 


where  q(g)   is  a  microfield  and  Z|^  is  the  partition  function  for  this 
system.     By  separating  out  the  integral  over  perturber  number  one,  we 
obtain 


t 


<f>     =— N_/d3r,f   r/d3(N-^V  6((Mj-  I     l.)e-^^\.  U.B.k) 

^       Z„Q(g)  '    '  '      i-2     ^ 

N 

Within  the  square  brackets,   V  can  be  considered  to  be  an  N-1   particle 
potential.     Then  [   ]  has  the  form  of  a  microfield  for  this   potential, 
times   its  partition  function, Z'    ,.     In  symbols,   [    ]  =  Q'(e-E,)Z'    ,. 
Note  that 


^=^  (I.B.5) 

^N  " 


where  g(r)  is  the  radial   distribution  function  and  n  is  the  system 

volume.     Now  <f>     can  be  written 
e 


<f>^   =  njd^r  f^g(r)Q'(E-E^)/Q(E)  (I.B.6) 

where  nj^   is  the  ion  density.     The  Independent  Perturber  Model    (Iglesias 
et  al.    1983)   is   used  for  Q'    and  Q.     See  Appendix  D.     In  this  model 


1       r  "^1^^) 

Q(e)   =  -^  J     dk  k  sin(ke)  e 

2tt  e     0 


where 


h   (k)   -  iJiT  J  dr  r^  e'^^(j    (kE(r))-1).  (I.B.7) 

0  " 


Therefore, 


Zen.  »  p  nh. (k)     » 

<f>     -  — ^ J     dk  k  e  J      dr  r  g(r)F(r)I     ,  (I.E. 8) 

^        2Tr'^Q(e)      0  0  " 


^.^ 


where  the  integral  over  angles,   I   ,   is  given  by 


Ijj  =  /dn(1-3u   )Jo(k|E-eJ).  (I.B.9) 


This  integral   is  calculated  in  Appendix  E  with  the  result 


I^  =  ^8ir  i^(ke)i^ikt^).  (I.B.10) 


Substituting  this  into  equation   (I.B.8)   yields 


00  f|h    (k )  '"' 

-STTZen.eJ  dk  k^e     ^       jpCke)/  dr  r^g(r)F(r)j    (ke    (r)) 

<f> i-° ,-j-r ° . 

J      dk  k  e  sin(ke) 

0 

(I.e. 11) 


This  can  be  expressed  in  terms  of  dimensionless  variables  as 


w    n  ^  -  -ax  2  2 


-e.j\,.h"^'^^2^'^'f^'-^'-^ 


<f>    = 

e 


^        (1    +  ax  +  — 2~^J2^-^^1^ 
0  " 

"nhTTkT 


J     dk  k  e  sin(ke) 

0 


(I.E. 12) 


where 
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2   4   3 
e  is  scaled  in  terms  of  Eq  =  e/r^  ,  -^^q   "g  *  "I » 

k  is  scaled  in  terms  of  Eq  , 

3 
f  is  scaled  in  terms  of  e/r.  , 

a  =  Tq/A  and 


Numerical   procedures  are  given  in  Appendix  F.     We  used  the  Debye-Huckel 
radial   distribution  function,   g(r),   in  this  calculation.     In  order  to 
estimate  the  error  introduced  in  this  way,  we  ran  one  case   (ng  = 
10^^  cm"^,   T  =  800  eV)   using  a  Monte  Carlo  generated  radial   distribution 
function.     At  the  value  of  c  corresponding  to  the  peak  of   the 
microfield,   the  difference  in  <^^^>^  was  less  than  5%. 

C.     Trends   and  Results 
As  seen  at  the  conclusion  of  Appendix  B,   the  quadrupole  interaction 
gives  rise  to  an  asymmetric  effect,  while  the  dipole  interaction  is 
symmetric.     A  measure  of  the  relative  importance  of  the  quadrupole 
interaction  on  the  line  broadening  problem  is  given  by  the  ratio  of  the 

quadrupole  interaction  to  the  dipole  interaction.     This  is  the  smallness 

2  2     1/3     4/3 

parameter  mentioned  earlier,      (n  aQ/Z)/R^Q  -■  n  ne       /Z 

That  is,  the  asymmetry  due  to  the  ion-quadrupole  interaction  will 

increase  with  principal   quantum  number  and  with  density  but  will 

decrease  with  radiator  nuclear  charge.     These  trends  can  be  seen  in  the 

calculated  profiles.     Figures   2a,   2b,  and  2c  show  the  Ly-a  lines  of 

hydrogenic  Argon,  with  and  without  the  quadrupole  interaction,  at 

densities  of   10^^,    3  x   10^^,   and  lO^'^  electrons  per  cubic  centimeter  and 


A 


V    ■;  V  ...  1 


FIGURE  2 


Ly-a  lines  of  kr*^'^ ,  with  and  without  the  ion-quadrupole 
interaction,  at   a  temperature  of   800  eV.     a)     Elec^^on  ^^ 
density  =  10^^  cm"'^;    b)     Electron  density  =3x10       cm 
c)     Electron  density  »  10       cra~-^. 
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a  temperature  of  800  eV.  These  figures  show  a  blue  asymmetry  (as 
predicted  in  section  I. A)  which  Increases  with  density.  Figures  3  and  H 
show  the  Ly-6  and  Ly-Y  lines  at  the  same  densities.  Again  the  increase 
in  asymmetry  with  density  is  observed.  By  comparing  Figures  2-ii   one  can 
also  notice  the  increase  in  asymmetry  with  increasing  initial  principal 
quantum  number.  Figure  5  shows  the  Ly-g  line  of  hydrogenic  Neon  at  an 
electron  density  of  3  x  10^^.   Comparison  with  Figure  3b  demonstrates 
the  expected  decrease  in  asymmetry  with  increasing  radiator  nuclear 
charge . 

The  quadrupole  interaction  in  line  broadening  is  the  major  thrust 
of  this  work,  however,  before  comparison  with  experiment  can  be  made, 
other  sources  of  asymmetry  must  be  investigated. 
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FIGURE  3 


Ly-B  lines  of  Ar"^^"^,  with  and  without  the  i o n- qua dru pole 
interaction,  at  a  temperature  of  800  eV.  a)  Electcon 
density  =  10^^  cm''^;  b)  Electron  density  =  3  x  10   cm 
c)  Electron  density  =  10   cm"^. 
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CHAPTER  II 
OTHER  ASYMMETRY  EFFECTS 

A.  Fine  Struct 'jre 

In  this  chapter  the  effects  of  fine  structure  and  quadratic  Stark 
effect,  as   presented  by  Woltz  and  Hooper   (1985),   are  discussed.     As  a 
starting  point,  we  take  equation   (I. A. 38).     The   quadrupole  interaction  may 
be  included  or  not  by  retaining  or  omitting  the  Q^^  term.     In  this  chapter 
we  will   be   interested     in  the  characteristic  features     of  the  other 
asymmetry  effects  and  so  will  not  include  the   quadrupole  interaction.     In 
Chapter  III,   all   three  effects  will   be   considered  together. 

First  we  consider  the  'ine  structure  effects.     The  Hamiltonian,  H^, 
of  equation   (I. A. 38)   is  often  taken  to  be 


H     =   -fi^V^/2M  -  Ze^/r  (II.A.1) 

r 


where  p  is  the  reduced  mass  of  the  electron  and  nucleus.     If  the 
Schrodlnger  equation  is  solved  with  this  Hamiltonian  the  eigenstates  are 
the  usual    |n!lm>  states  where  n  is  the  principal   quantum  number,  I  is  the 
orbital   quantum  number,  and  m  is  the  magnetic  quantum  number.     The 
eigenvalues,   given  by 


<nam|H    |nam>  -  -  Hl^  ^  (II. A. 2) 


are,   for  a  given  principal   quantum  number,  all  degerate.     Consider  the 
states  with  a  given  principal   quantum  number.     If  a  uniform  electric  field 

-34-  .  . 
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is  applied,  the  new  energy  levels  can  be  derived  from  the  linear  Stark 
effect,  which  is  the  result  of  first  order  degenerate  perturbation 
theory.     Some  of  the  states  are  raised  in  energy  and  others  are  lowered, 
maintaining  the  symmetry  about  the  original  energy  level.     This  symmetry 
manifests  itself  in  a  symmetric  line  profile.     If  the  original  Hamiltonlan 
did  not   give  rise  to  completely  degenerate  states,  an  asymmetry  would  be 
introduced.     The  degeneracy  is  actually  lifted  by  two  effects  neglected  in 
equation   (II.A.1).     the  relativistic  effects  and  the  spin-orbit 
interaction.     These  effects  are  of  the  same  order  of  magnitude  and 
together  are   known  as  fine  structure.     The  exact  solution  of  the  Dirac 
formula  gives  a  result  which  is  diagonal   in  the  representation   |nljj^>. 
The  result   is    (Bethe  and  Salpeter    1977), 


<nUjjHjn^jJ^>  =   [1    ^   [^ ^     3^-   ^'^"''''  '   '• 

n-(j+1/2)  +  /(j+1/2)  -a   Z         (II. A. 3) 

where  a  is  the  fine  structure  constant  (-1/137).  and  j  is  the  quantum 


num 


ber  associated  with  the  total  angular  momentum.     Little  error  is 


2 

introduced  in  expanding  this  result  to  first  order  in   (aZ)    , 


<n«J  jH^|nWJ,>  -  ^  ^[1  *  ir-(FT75  '  i^'l-      '"•'•"' 


The  first  term  is  Just  the  result  given  in  equation   (II. A. 2)  and  the 
second  term  is  the  fine  structure  correction.     The  Dirac  theory  does  not 
contain  radiative  corrections,  but   these  are  of  order  aln(a)  smaller  than 
the  last  term  kept  above. 


If  we  calculate  JCu.e)   using  the  basis   |nilni  m  >,  we  have 


Is 


1 

J(a),e)   = Im     I  <nilm  m    |D|n£'m   'ra   '> 

I     s 


0  <nj,'ni   'm   '  |[a)-H  -a)-.-p'e  -5'/((d)]      |nilm  m  > 


=  -I  ImTr  D^-'^R^^^  (II. A. 5) 

IT 

r 


If  we  neglect   nonradiati ve  transtitions   between  states  of  different 
principal   quantum  number  and  consider  only  the  linear  Stark  effect,  the 
matrix  R^"^^  is  obtained  by  inverting  the  matrix  having  elements 


<n)lm„m    L-H  -to  -p«e-  ^(oo)  Ini'm. 'm  '>.  (II. A. 6) 

S,  s '       r     0  '  S,     3 


The  operator  m-m-p-t-^iui)   is  usually  calculated  in  an   |nilm  > 
representation.     The  only  difference  here  is  that   it  will   have  twice  as 
many  rows  and  columns    (for  the  two  values  of  fflg),  with  the  elements   being 
equal  to  zero  if  m     *  m   ' .       Equation   (II. A. 3)     gives  H_  in  the 
representation      |nljj   >.       It  is  a  simple  matter  to  transform  to  the 
Inlm.m  >  basis  using  the  unitary  treuisformation 

At       S 


-11+1/2-J  i        1/2  j 


(II. A. 7) 


K. 


!p/*-J'l"-  '  <S^.  •  V  -   V.^.r; 


The  symbol 

is  the  Wigner   3-j   symbol  and  is  easily  calculated   (Edmonds   1957). 
We  are  now  able  to  calculate  lines  with  fine  structure.     For 
convenience,  we  take  the  energy  level  of  the  state  with  the  largest  value 
of  j   to  be  our   zero  reference  point.     Figures   6a,     6b,  and  5c  show 
the  Ly-a  lines  of  hydrogen! c  Argon,  with  and  without  fine  structure,   at 
densities  of   10     ,    3  x   io     ,   and   10^     electrons   per  cubic  centimeter  and  a 
temperature  of  800  eV.     Figures   7  and  8  give  the  same  cases   for  the  Ly-B 

and     Ly-T  lines.     Figure  9  shows  a  Ly~B  line  of  Neon  at   a  density  of 

23 
3x10       electrons  per  cubic  centimeter. 

From  equation   (II. A. 4),  we  see  that  the  fine  structure  correction, 

and  therefore  the  asymmetry  of  the  line  shape,   increases  with  nuclear 

charge,   Z,   but   decreases  with  increasing  principal   quantum  number,  n. 

Comparison  of  Figure   9  and  Figure   8b  shows  the  expected  dependence  on 

nuclear   charge  and  a  comparison  of  the  appropriate  members  of  Figures   5-8 

shows  decreasing     asymmetry  with  increasing  principal   quantum  number.     The 

decrease  in  asymmetry  with  increasing  density   (see  Figure  5,    7,  or   8)   is 

also  to  be  expected.     As  the  density  increases,  so  does  the  average  ion 

field  perturbation  of  the  energy  levels.     The  fine  structure  correction  is 

purely  a  result  of  atomic  Physics  and  does  not   change  with  density.     At 

low  densities,  the  fine  structure  correction  is  larger  than  the  average 

ion  field  perturbation  and  therefore  a  large  asymmetry  exists.     At  high 

densities  the  reverse  is  true  and  the  line  is  more  symmetric. 


y 


FIGURE   6 

Ly-a  lines  of  Ar''^'^,  with  and  without  fine  struct ^e,   at  a 
temperature  of   800  eV.     a)     Electron  density  =10       cm     ; 
b)     Electron  density  =   3  x   lO''^  cm  ^   ;   c)     Electron  density 
=  10^^  cm'^. 
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FIGURE   7 

Ly-6  lines  of  Ar*^"^,  with  and  without  fine  structure  at  a 
temperature  of   800  eV.     a)     Electron  density  =  10  ^  cm~^; 
b)     Electron  density  =   3  «   1 0  ^  cm  ^   ;   c)     Electron  density 
=  10^^  cm"3. 
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FIGURE  8 


+  1  7 

Ly-Y  lines  of  Ar   ,  with  and  without  fine  structure  at  a 
temperature  of  800  eV.  a)  Electron  density  =  10-^  cm~-^; 
b)  Electron  density  =  3  x  1 0   cm   ;  c)  Electron  density 
=  lO'^^  cm  ^. 
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B.     Quadratic  Stark  Effect 
The  quadratic  Stark  effect   causes  asymmetry  in  hydrogenic  spectra. 

This  effect  is  one  order  higher  in  our  smallness  parameter, 

2 
(n  a^/D/R     ,   than  the   quadrupole  interaction   (Sholin   1969). 

Consider  again  the  Hamiltonian   (II.A.1)     for  an  isolated  hydrogenic 

ion  of  nuclear  charge  Z.     This  Hamiltonian  has  eigenstates    |nlm>  and 

eigenvalues 


'nam       __2       2 
2fi       n 


(II. B.I) 


which  are  completely  degenerate  in  i   and  m  .  If  a  uniform  electric  field 
in  the  z  direction  is  applied,  the  total  Hamiltonian  is  now  H  =  Hp  +  H^ , 
where  H^  =  eze,  e  being  the  magnitude  of  the  electric  field.  Degenerate 
perturbation  theory  provides  a  systematic  procedure  for  finding  the  new 
energy  levels.  To  first  order,  this  procedure  amounts  to  diagonalizing 
H^  ,  within  each  block  of  degenerate  states,  that  is,  within  each  block  of 
states  having  a  given  principal  quantum  number.  The  results  are  the  usual 
Stark  states  |nqm>  with  energy  levels  given  by 


2n  n 


Higher  order  corrections  are  increasingly  more  difficult  to  calculate. 
They  involve  the  mixing  of  states  with  different  principal  quantum  numbers 
and  therefore  require  infinite  sums.  The  sums  for  the  second  order 
correction,  the  quadratic  Stark  effect,  can  be  done  analytically  (Bethe 
and  Salpeter  1977)  with  the  result 


m 
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^nli   =  ^nl!  -  T6  a,3(|)'^e2(17n2-3q'-9m2.i9)        (II.B.3) 
nqm    ncjm   i  b  0  Z 


where  Hq  is  the  Bohr  radius. 

Consider  now,  J((d,c)  as  given  by  equation  (I. A. 38),  neglecting,  for 
the  moment,  the  quadrupole  term.  Let  us  refer  to  the  operator  that  is 
to  be  inverted  as  the  resolvent  matrix.  If  we  consider  e  constant, 
J((i),e)  as  a  function  of  co  has  the  form  of  a  Lorentzian,  whose  peaks 
correspond  to  the  hydrogenic  energy  levels  for  that  fixed  e.   A  common 
procedure  is  to  numerically  invert  the  resolvent  within  the  subspace  of 
initial  radiator  states  (those  states  with  principal  quantum  number  = 
n).  This  is  clearly  equivalent  to  the  linear  Stark  effect  of  degenerate 
perturbation  theory,  which  was  outlined  above.  For  this  reason,  the 
approximation  of  tracing  only  over  initial  radiator  states  is  known  as  the 
linear  Stark  effect.  This  will  be  a  good  approximation  as  long  as  the  ion 
field  perturbations  are  small  compared  to  the  separation  between  states  of 
different  principal  quantum  number. 

It  is  not  clear  how  to  systematically  improve  the  line  broadening 
theory  to  include  the  quadratic  Stark  effect.  We  shall  outline  here  two 
quite  different  approximate  treatments  of  the  quadratic  Stark  effect  in 
line  broadening.  These  two  different  approximations  give  results  that  are 
in  good  agreement,  suggesting  that  both  work  well  for  the  plasma 
conditions  discussed  here. 

The  states  that  are  closest  in  energy  to  the  states  with  principal 
quantum  number  n,  (and  therefore  most  likely  to  mix  with  them)  are  those 
states  with  principal  quantum  number  n+1 .  The  first  method,  then,  is  to 
include  in  the  calculation  of  the  resolvent,  not  only  all  states  with 
principal  quantum  number  n,  but  also  n+1.  The  second  method  uses  Stark 


r■f■..^' 


states  |nqm>  of  only  one  principal  quantum  number,  but  artificially  uses 
the  quadratic  energy  shift  (II. B. 3)  instead  of  the  linear  term  (II. B. 2) 
for  the  diagonal  matrix  elements  of  Hp  +  H^ . 

Figure  10a,  b,  and  c  show  Ly-a  lines  of  hydrogen! c  Argon,  with  and 

without  the  quadratic  Stark  effect,  at  densities  of  ^0^^,    3  x  10^^,  and 

24 

10  electrons  per  cubic  centimeter  and  a  temperature  of   800  eV.     Figures 

11  and   12  give  the  Ly-g  and  Ly-Y  lines  for  the  same  conditions.     Figure   13 
shows  the  Ly-6  line  of  hydrogenic  Neon  at   3  x   10       electrons  per  cubic 
centimeter.     From  equation   (II. B. 3),  we  see  the  effect  of  the  quadratic 
Stark  effect  should  increase  with  principal   quantum  number  and  with 
density   (since  increasing  the  density  increases  the  average  ion  field)  and 
should  decrease  with  increasing  nuclear  charge.     These  trends  are  born  out 
by  the  figures. 

The  Lyman-B  line  is  particularly  Interesting  because  it  is  usually 
intense  enough  to  be  seen  when  the  Lyman-a  line  is  observable,  but   it  is 
not  as  optically  thick.     The  asymmetry  of  the  Lyman-6   induced  by  the 
quadratic  Stark  effect   can  be  explained  as  follows:     In  the  linear  Stark 
effect,  the  electric  field  causes  the  initial  states  of  the  Lyman-B  line 
to  be  split  symmetrically  above  and  below  their  unperturbed  energy  level, 
the  splitting  being  proportional  to  the  field  strength.     The  quadratic 
Stark  effect  slightly  lowers  the  energies  of  all  these  states;   and  the 
lowering  becomes  greater  as  the  field  strength  increases.     This  gives  rise 
to  an  increased  density  of  states  contributing  to  the  blue  peak  of  the 
line  profile  and  a  decreased  density  of  states  contributing  to  the  red 
peak,  causing  a  blue  asymmetry. 

It  is  a  simple  matter  to  combine  the     quadratic  Stark  effect,  as 
outlined  here,  with  the  fine  structure.     In  figure   14,  we  show 


FIGURE  10 

Ly-a  of  Ar*^"^,  with  and  without  the  quadratic  Stark  effect,  at 

a  temperature  of  800  eV.  a)  Electron  density  =  10  ^  cm'^; 

b)  Electron  density  »  3  x  1 0  -^  cm~  ;  c)  Electron  density 
-  10^^  cm'3. 
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FIGURE   11 


Ly-6of  Ar"^^^,  with  and  without  the   quadratic  Stark  effect,   at 
a  temperature  of   800  eV.     a)     Electron  density  =  10  ^  cm~-^; 
b)     Electron  density  =   3  x   10       cm'     ;   c)     Electron  density 
«  10^^  cm~3. 
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FIGURE    12 
+  1  7 

Ly^T  of  Ar        ,  with  and  without  the   quadratic  Stark  effect,  at 
a  temp)erature  of   800  eV.     a)     Election  density  =  10-^  cm"-^; 
b)     Electron  density  =3x10       cm       ;    c)     Electron  density 
=102^cm-3. 
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FIGURE    14 
JCoj.g)  of  a  Ly-ci  line  for   various   values  of  £, 
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J((ij,e),  for  a  Ly-a  line,  as  a  function  of  to,  for  various  values  of  ion 
field,  e.   The  peaks  should  correspond  to  the  energy  levels  of  the 
hydrogenic  states  for  different  values.  In  Figure  15  we  show  the  energy 
levels  of  the  n  =  2  states  of  hydrogen  as  a  function  of  electric  field 
(Luders  1951).  Comparison  of  Figures  14  and  15  show  good  agreement, 
giving  confidence  in  the  methods  outlined  here. 


^fii 


FIGURE   15 

Energy  levels  of  the  n=2  states  of  hydrogen  as  a  function  of 
perturbing  field  strength. 
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CHAPTER  III 
RESULTS 

It  is  a  simple  matter  to  combine  the  three  asymmetry  effects  dis- 
cussed in  the   previous  chapters.     The   quadrupole  interaction  is  included 
in  J(a),e)  as   in  equation   (I. A. 38).     Fine  structure  is  included  by 
calculating  Hp  as   in  equation   (II. A. 4).     The  resolvent  matrix  is 
calculated  only  for  states  of  a  given  principal   quantum  number.     It  has 
2n     rows  and  columns,  twice  the  number  it  would  have  if  spin  were 
ignored.     This  result   can   be  transformed  to  the    |nqm>  representation  so 

that  the   quadratic  Stark  effect  may  be   included.     In  this  representation, 

■*■  -*■ 
the  operator  p'e  is  assumed  diagonal  and  is  replaced  by 


|^-i^a„3,|,V<UnW-9.^-19)  (ni.».„ 


as  in  equations    (II. B. 2)  and   (II. B. 3). 

In  Figures   16-18  we  show  line  shapes  with  and  without  all  three 
asymmetry  effects  included.     Figure   1 6a ,  b,  and  c  show  Ly-a  lines  at 
densities  of   10^^,    3x10"^,   and   10^     electrons   per  cubic  centimeter  and  a 
temperature  of  800  eV.     Figures    17  and   18  show  the  Ly-6  and  Ly-Y  lines 
under  the  same  conditions. 

As  mentioned  earlier,  Lyman-6  lines  are  particularly  useful  as 

density  diagnostics   because  they  are  nearly  as   intense  as  Lyman-a  lines, 

■  but  are  not  nearly  so  optically  thick.     Although  inclusion  of  the  three 

asymmetry  effects   discussed  here  will  cause  significant  changes  in  Lyman 

line  profiles,  for  most  temperature  and  density  regimes   the  Lyraan^B  line 

'%:'-■  .  ..  '  -62-  '      '       ■ 


FIGURE  16 


Ly-a  lines  of  Ar"^^"^,  with  and  without  the  three  asymmetry 
effects,  at  a  temperature  of  800  eV.  a)  Elect ron^,   ,3 
density  =  10^-^  cm~3;  b)  Electron  density  =  3  x  1 0   cm 
c)  Electron  density  =  10   cm~3. 
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FIGURE  17 


Ly-B  lines  of  Ar"^^"^,  with  and  without  the  three  asymmetry 

effects,  at  a  temperature  of  800  eV.  a)  Electron^^  _^ 

density  =  10^^  cm"^;  b)  Electron  density  =  3  x  10   cm 
c)  Electron  density  =  10   cm~^. 
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FIGURE  18 


Ly-Y  lines  of  Ar"*"^"^,  with  and  without  the  three  asymmetry 
effects,  at  a  temperature  of  800  eV.  a)  Electron     ^ 
density  =  10^^  cm  ^;  b)  Electron  density  =«  3  x  i  Q   cm 
c)  Electron  density  =  10   cm~-^. 
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will  still  be  a  two  peaked  function.  Let  us  define  a  measure  of  asymmetry 
for  the  Lyman-B  line  as  2(1^^  -Ip)/(Ifj  +1^)  where  I^^  and  Ip  are  the  blue 
and  red  peak  intensities,  respectively.  In  Figure  19,  we  plot  this 
asymmetry  as  a  function  of  density  for  an  Ar*  "^  plasma  at  800  eV.  Lines 
which  include  only  the  fine  structure  splitting  have  a  negative  (red) 
asymmetry  which  decreases  in  magnitude  with  increasing  density.  Both  the 
quadratic  Stark  effect  and  the  ion  quadrupole  interaction  cause  positive 
(blue)  asymmetries  which  increase  with  increasing  density.  These  two 
effects  are  shown  together  in  Figure  19.  Also  shown  is  the  asymmetry  of 
lines  which  are  calculated  using  all  three  effects.  The  asymmetry  of 
these  lines  is  a  rapidly  changing  monotonic  function  of  density  (at  least 
in  the  range  of  experimental  interest)  and  therefore  could  provide 
diagnostics  if  the  effects  of  opacity  are  not  to  great.  It  is  interesting 
to  note  that,  while  each  effect,  taken  alone,  would  cause  significant 
asymmetry  at  4  x  IQ  -"  cm  ,  their  combined  effect  is  to  yield  a  nearly 
symmetric  line  profile. 

Since  this  is  the  region  of  current  experiments,  these  asymmetries 
will  be  very  difficult  to  observe.  The  table  below  shows  the  separate 
contributions  to  the  asymmetry  at  several  densities. 
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TABLE  1.  Contributions  to  the  asymmetry  of  the  Ly-6  line  of  Ar"^^*^ 


10 


23 


3x10 


23 


10 


24 


3x10 


24 


FINE  STRUCTURE 

-.187 

QUADRATIC  STARK 

.023 

QUADRUPOLE 

.028 

TOTAL 

-.142 

.112 

-.053 

-.030 

.034 

.062 

.064 

.051 

.098 

.151 

.022 

.085 

.174 

If  this  procedure  is  performed  for  different  elements,  the  trends  are 
the  same,  but  the  point  of  zero  asymmetry  will  be  different.  For  example, 
Ne"^^  lines  are  nearly  symmetric  at  10^^  cm"^. 

The  asymmetry  of  experimental  emission  line  profiles  can  be  affected 
by  more  than  just  the  three  effects  outlined  here.  For  example,  the 
background  radiation  slope  must  be  suitably  accounted  for  and  pertinent 
overlapping  spectral  lines  must  be  included  in  the  calculation  (Delamater 
1984). 


CHAPTER   IV 
SPECTRAL  LINES  FROM  HELIUM-LIKE  IONS 

It  has   been  conjectared   (H.R.Griem,  private  communication)  that  for 
high  temperatures  and  densities,   a  product  of  hydrogenic  wavef unctions 
would  provide  a  satisfactory  approximation  for  the  helium-like  wave- 
functions  of  high-Z  ions.     A  strict  interpretation  of  this  statement, 
however,  would  imply  a  degeneracy  in  the  orbital   quantum  number,    "I",  a 
degeneracy  which  is  not  present  is  helium-like  ions.     This  non-degeneracy 
causes  an  asymmetry  in  the  helium-like  lines  even  without  considering  the 
ion  quadrupole  interaction,   the  fine  structure  splitting,   or  the   quadratic 
Stark  effect  which  give  rise  to  additional   asymmetries,   as  seen  in  earlier 
chapters.     Because  of  this  qualitative  difference  between  the  hydrogenic 
and  helium-like  wavef unctions ,  we  take  the  following  as  our  model.     One 
electron  is  to  be  in  the  n=1   level  and  the  other  initially  is  in  a  higher 
level.     The  ground  state  electron  has  only  two  effects    (in  our 
approximation).     First,   it  causes   the  i-degeneracy  for  the  other  electron 
to  be  broken,  and,   secondly,  and  less  importantly,   it  screens  the  nucleus 
so  that  matrix  elements  are  calculated  using  hydrogenic  wavef unctions  with 
a  nuclear  charge  of  Z-1 . 

Initially,  the  two  electrons  may  couple  into  either  a  singlet  or  a 
triplet  state,   that  is,   a  state  of  total  spin  zero  or  one.     The  ground 
state,  however,  must  be  a  singlet  by  the  Pauli  exclusion  principal. 

Transitions  from  triplet  states  to  singlet  states,  such  as 

3  2  1 

1s2p     P  -»•   Is       S,     give  rise  to  lines  called  intercombination  lines. 

These  lines  have  low  transition  probabilities  since  they  are  dipole 
:»^•^  ,  -73- 
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forbidden.     For  this  reason,  we  will   completely  ignore  the  triplet  states 
in  this  calculation  and  only  consider  singlet-singlet  transitions.     In 
actual  experiments,   the  intercombination  lines  may  be  observable.     In 
these  cases  they  may  be  approximated  by  a  Lorentzian  or  Voight  profile 
(Delamater   1984). 

Following  the  formalism  of  Chapter  I,   but  neglecting  the  quadrupole 
interaction,  we  have 

I((d)    -  Jde  Q(e)J(u,e), 

J(a),e)    =  -  ;J-  Im  Tr  dU-CH^-Uq)   -  p«£  -  P^  U)}~\  (IV. A.I) 

a 

The  trace  is  taken  with  respect  to  singlet  states  of  a  given  principal 
quantum  number.     The  Hamiltonian,   Hj^.is  assumed  to  have  eigenvalues 

2 

^n&m   °  "   ^^~2^      *  ai(n.il,Z).  (IV. A. 2) 

n 

where  the  energy  has  been  expressed  in  Rydbergs .     The  first  term  is  just 
the  energy  of  an  isolated  hydrogenic  ion  of  nuclear  charge  Z-1   and  the 
second  term  is  an  i-dependent  energy  shift.     This  shift  represents  the 
effect  that  the  ground  state  electron  has  on  the  energy  levels  of   the 
excited  electron.     Clearly,   u)(n,X,,Z)  removes  the  ^-degeneracy. 

The  values  for  a)(n,£,Z)  must  be  input  to  our  code  from  an  outside 
source.     For  the  helium-like  a-lines,   it  is  the  n«2  helium-like  energy 
levels  that  are  required.     These  have  been  conveniently  calculated  by 
Knight  and  Scherr    (1963)   as  an  expansion  in  powers  of   (1/Z).     For  higher 


Ite-:- 
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series  members,  the  energy  levels  must  be  obtained  elsewhere,  either  from 
theory  or  experiment   (Bashkin  and  Stoner   1975;   A.Hauer,   private 
commun  i  cat  i  on ) . 

If  the  ion  field  perturbation  is  small  compared  to  the  energy 
difference  between  states  of  different  orbital   quantum 
number,   a)(n,il^,Z)  -  a)(n,)l2.Z).   then  these  states  will  not  mix,   and  there 
will   be  no  linear  Stark  effect.     The  energy  levels  will  shift 
quadratically  with  electric  field.     For  large  ion  field  perfurbations ,   the 
states  mix  substantially  and  the  levels  shift  linearly  with  electric 
field.      Appendix  G  shows  this  change  from  quadratic  to  linear  dependence 
for  the  simplest  case,   the 
a-line. 

It  is  interesting  to  note  that,   in  the  case  of  the     a-line  given  in 

Appendix  G,   the  shifted  components  shift  symmetrically  about  their  initial 

mean,   -u     /2,   and  yet  the  line  is  markedly  asymmetric.     There  are  two 
sp 

reasons  for  this  asymmetry.     First,   the  unshifted  components  are  centered 

about   zero   (the  i!,=1   energy  level)  rather  than  the  symmetry  point 

-oj     /2.       Secondly,   the  transition  probabilities  for  the  shifted 
sp 

components  are  not  equal.  The  shifted  components  are  linear  combinations 
of  the  states  |200>  and  |210>.  In  the  absence  of  any  field  these  states 
are  not  mixed  and  have  relative  energy  levels  -ui  and  zero.  This  gives 
rise  to  asyramerty  because  the  state  |200>  is  dipole  forbidden.  Only  for 
large  electric  fields  is  the  situation  symmetric,  as  the  eigenvalues  are 
then   (|200>  ±    |210>)//2. 

While  the  presence  of  more  different  initial  energy  levels  introduces 
additional   asymmetry  to  the  lines  of  higher  series  members,   it  is  the  a 
line  which  is  qualitatively  the  most  different  from  its  hydrogenic 
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counterpart.     Figures   20-22  show  hellom-like  a,   6,  and  Y  lines  of  Argon  in 
regimes  of  experimental   interest.     Because  the  a  line  is  so  narrow,  we 
show  it  on  a  log  scale  to  bring  out  the  features. 

The  ability  to  calculate  helium-like  spectra,   as  well  as  hydrogenic 
line  shapes,  greatly  increases  the  possibility  of  using  several  lines  to 
analyze  a  single  experimental  shot.     A  problem  with  this,  however,   is  that 
the  peak  hydrogenic  line  intensity  and  the  corresponding  peak  helium-like 
line  intensity  will,   in  general,  occur  at  different  temperatures  and 
densities  that  arise  during  the  course  of   the  experiment.     Therefore  the 
choice  of  temperature  and  density  conditions  which  best  fit  a  hydrogenic 
line  in  a  given  experiment  will  not  necessarily  fit  a  helium-like  line 
occuring  in  the  same  experiment.     This   problem  may  be  greatly  reduced 
through  the  time-resolved  spectroscopy  techniques  presently  available. 


FIGURE   20 


Helium-like  a-line  of  Ar"^    '   at  an  electron  density  of 
10^  cm~-^  and  a  temperature  o*"   800  eV. 
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FIGURE   21 


Helium-like  B-line  of  Ar"^^"^  at  an  electron  density  of 
10^  cm"-^  and  a  temperature  of   800  eV. 
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FIGURE   22 

Helium-like  Y-line 
10^  cra'^  and  a  tei 


of  Ar"^^^  at  an  electron  density  of 
emperature  of  800  eV. 
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APPENDIX  A 
ZWANZIG  PROJECTION  OPERATOR  TECHNIQUES 


From  Eqn.  (I.A.I  5)  we  have 


Ko))  =  JdH'  W(4')J(a),4'),  (A.l) 


and 


JCiu.'i')  =«  tt'  Re  J  dt  e^  Tr  d'f  (a  ,f  )<d(-t  )>, 
0       a 


where  <X>  =  f"  (a,y)  Tr  p('i), 


'(a,4')  =  Tr  p  6(4'-$  ), 


W(«i')  =  Tr  p  (5(4'-*  ), 


and 


p(4')  =  p  5(4'-*  )/W(4'). 


Define  an  operator  P  by 


PX  =   <X>  =  f''\a,4')  Tr  p(f)X.  (A. 2) 


Note  that 
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.2,.    -1,_  ....  _  ^.,...^.-1 


P^X  -  f"  (a.f)  Tr  [p(4')f~  (a.H-)  Tr  p(T)x] 


=  f^''(a,l')(Tr  pCf))  f'^Ca.t)  Tr  p(4')X 


PX  (A.3) 


so  P^  =  P. 


This  property,  idempotence ,  defines  a  projection  operator  (Reed  and  Simon 
1972).  We  can  now  write  J(a),'l')  as 


JCoj.^-)  =  Ti'^Re  Tr  d-f  (a,4')  j  dt  e^*^  Pd(-t) 
a  0 


ir"''Re  Tr  d-fCa.f)  p5^(co),  (A.'O 


where  c^,  (oj)   =  pA(w)  ,   and  ^(lo)  =   /  dt  e  '^  d(-t)     The  Zwanzig  Projetion 

^  0 

operator  technique  can  be   used  to  obtain  a  general  expression  for 

(?^^  (w) . 


Let  D(t)   =  d(-t).       We  have 


i  5(t)   =  LD(t).  (A. 5) 


Let  PD  =  D^  and  (l-P)D  =  D2  so  that  D^  +  D2  =  D.  Then 
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i  D^  =  (1-P)L(D^  +  D^), 

and  ( A . 6 ) 

i  D^  =  PL(D^  +  D^). 


Therefore 


(3^  +  i  (1^P)L)D2(t)  =  ^i(l-^P)LD^(t),  (A. 7) 


which  has  the  solution 


D  (t)  =  e^^^^^'^\    (0)   -  i  J  e'^^^^"^^^(1-P)LD  (t-s)d3.   (A. 8) 

0 


For  ovJT  choice  of  P,   D2(0)   =  0.     Therefore 


i  5    (t)   =  PLD    (t)  -   i   /     PL  e''^^^^**''^^(1-P)LD    (t-s)d3.  (A. 9) 

'  '  0 


This  can  be  transformed   using  the  Convolution  Theorem, 


,    r       iujt*    ,....        „,     f     ia)t_    ,....        ,„,    c       iuot  -'it(1-P)L 
i   J      e       D.(t)dt   =   PL  J    e       D,(t)dt  -   iPL  J      e       e 


dt 


or 


9    (1*P)L  J     e^'^^D    (s)ds,  (A. 10) 

0 


-iD,(0)    +  i^^Aoi)   =   PlJ&.Coj)    +  PL  — ,]    -.,  (l-P)Lo^,  (a>).  (A. 11) 

1  1  1  a)-U-P)L  1 


So 

(0^(0))   =   i[u)-PL'-PL((i)^QL)'^QL]''^D^(0),  (A. 12) 


!/■ 
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where  Q  =   1-P.     Putting  this  in   (A. 4)   yields 

J(a),4')   =  -ir*"   Im  Tr   d^f  (a.f)    [a)-<L>   -   <L(a)-QL)"S]d.  (A. 13) 


Let  L  =  L.  +  L^C*  )  and  L^  =  L(a,4')  +  L„  where  Ln  is  the  Liouville 


operator  for  the  isolated  atom  and  Lp  is  the  Liouville  operator  for  the 
plasma.  Then, 


<L>  =  <L(a,t(;)>  +  <L  >  +  <Lj.((f*)>.  (A.U) 


The  first  term  is 


<L(a,4')>X  =  f'^a,*)  Tr  p(f  )L(a  ,'1' )X  =  L(a,'}')X   , 
a  p  a        a 


and  the  second  term  is 


<L  >X     =  f'^a,*)   Tr  p(T)(L  X   )   =  0, 
p     d  pa 


because  LpX^  =  0,   so     <L>  =  L(a,'}')   +   <L^ii>   )>     when  acting  on  atomic 
functions.     Also 


P 


L^X^  -  f'^a.l-)   Tr  p('l')L^Xg 
P 


Va^^a^a  =  °-  (A- 15) 
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This  allows  us  to  replace  the  last  L  by     L  (*   )     in    ^(to).       Now, 
using  also  Q^  =  Q,  we  can  write 


» 


Now 


J(a),4')  -  -ir'  Im  Tr  d-f  (a,M')[u)-L(a,4')-<Lj(*  )> 

a 

-<L(unQLQ)~^QLj(**)>]~^  d. 


we  work  on  the  first  L  in  ^(w).   Consider 


PLQX  =  f~^  (a.y)  Tr  p('{')LQX 


f~\a,y)  Tr  Lp(4')QX 


(A.16) 


-■^(a,*)  Tr  (LCa.f)  +  L  +  Lj(*  ))p(4')QX.  (A. 17) 


The  first  term  is 


f'^Ca.f)  Tr  L(a,4')p('l')QX 


f'Ua,'9)Uaj)   Tr  p(4')X 


-  f  '^|a,4')L(a,4')  Tr  p(4')f  "la,f )  Tr  p('i')X 
P  P 

-  f"'^(a,«i')L(a,y)  Tr  p(*)X  -  f"""  (a,4')L(a,4')  Tr  p(Y)X  =  0.    (A. 18) 
P  P 

The  second  terra  on  the  right  hand  side  of    (A. 17)   is      ■    '  ■ 
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r'\a,^)   Tr  L  p(4')QX  «  f'Va.H-)   Tr  L  p(4')QX. 
p        P  p        P 


But     Tr  L„M  =  Tr  H  M  -  MH     .       Since  H^,  is  an  operator  only  on  perturber 
p       P         p       P  P  P  y         i^ 

states  and  these  states  are  being  traced  over,   it   is  legitimate  to 

cyclicly  permute  H_,,   i.e.     Tr  H  M  =  Tr  MH   .       Therefore, 

p  p  p 

P  P         ^ 

f'^a.f)   Tr   L   p(4')QX   =   0.  (A.  19) 

P 


This  leaves 


J(a),4')   =  -ir~     Im  Tr  d«f  (a  ,y  )[a)-L(a  ,4' )  -   <L   ($   )> 


_  .'1/_   ...s   „„   r    ,.»^    ,..,x,      ^.„^-U.     ,.».i-1 


(a,f)   Tr   Lj($    )p  (f)  (oi-QLQ)"  QLj($    )]"  d.  (A. 20) 

P 


This  is  the  result  shown  in  equation   (I. A. 18) 
Now  consider 


(a),'?)   =  f"^(a)   Tr  L^p  (¥)  (unQLQQ)~^QLj 
P 


f'"''(a)   Tr  r  pCV)  -     I      (-  QL-Q)"ql^.  (A. 21) 

p       ^  <"  n=0     "       °  ^ 


From  equation    (A. 18)   we  have  PLqQX  =  0,   or  QLqQX  =   LqQX  for   arbitrary  X, 
Using  this,  one  may  write 


OD 

(o).'}')   =  f'"\a)   Tr  L  pCf)  -     I     [I  i  qfqi 

I  w        „     0)     0 


to 
n=0 


""^a)   Tr  Ljp(4')(co-Lq)'^QLj.  (A. 22) 

P 


'^' 
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Noting  that  (a),¥)   acts  only  on  atomic  functions,  we  can  write 


f~    (a)   Tr  p.^Cy-**)    Tr  p     L,.X 
1  ae  I   a 

PLjX^ 5—^ =   0.  (A.23) 

Tr   p.  6  (4'-*   )  Tr   p 
I  ae 

This  yields 


(uj.f)   =  f'^a)   Tr  Ljp(4')(a)-LQ)~^Lj    .  (A. 24) 

P 


which  is  the  result  shown  in  eqn.    (I. A. 25). 


APPENDIX  B 

PROOF  THAT  THE  OCTAPOLE  EFFECT   IS  SYMMETRIC  TO 

FIRST  ORDER   IN   PERTURBATION 

In  this  appendix,   it  will   be  shown  that  terms  of  the  multipole 
expansion  with  odd  "X,"   values   contribute  a  symmetric  effect   on  the  line 
shape,  to  first  order  in  perturbation  theory. 

Consider  hydrogenic  states    |nqm>  in  the   parabolic  representation. 
The  presence  of  a  uniform  electric  field  removes   the  q  degeneracy  but   not 
the  m  degeneracy.     Upon  perturbation  by  an     J,m  mode,   the  energy  shift   of 
levels  with  given  n  and  q  values  will   be  a  linear  combination  of  terms 
like 


M  =  <nqm    |r  Y |nqm  >. 


8,m 


(B.I) 


Switching  to  the  spherical  representation  by 


|nqm>  =  I    |nj,m><n£m|ncp> 
X, 


where 


X   «    I        V        /    ,\1 /2(1+m*»q'*n) /TT-TT 
<nX,m|nqm>   =   ('-1)  ^       /2)!,+1 


f  nil  elL     0 

2  2  '^ 


m-q     m-t-q     ^ 


m 


(B.2) 


gives 
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M=     I       <ni^mJrX_|nJl2m2>(-1) 


1/2(1 +m^-q-n) 


^1^2 


1/2(1+M  -q-n)  

(-1)  /(2£^  +  1)(2il2+1) 


ml    ml  5. 

2  2  ^1 

m^-q     ra^+q 

—2 T  ""l, 


n-1        n-1 
2  2 


ra^-q     m^+q 

-2 r    ''"2; 


(B.3) 


Now, 


<ni  m    |r^_Jnl  m  >  =  J  R       R         r^*  dr  Jdfi  Y         Y__Y  . 

'    '  Im  0     "^1   ""^2  ^11    JLm  *2'"2 


(B.4) 


The  angular   integral   is   (Edmonds    1957) 


m,      (21   +1)(2il   +1)(2i+1)    /  ^1      ^     ^2]      /    ^1        ^        ^2  \ 


(-1)  V-::i 


4ir 


(B.5) 


0       0     0  7      \  -m,        m       my 


Let   R  represent  the  radial   integral. 


m^.m 


m.  (-^5-^)-q-n  (2)1+1)^(211+1)^(21+1) 

M  =     I     R  (-1)    \-^)        2  /     — ^ j^ 

Z^Jl^     12  ^"^ 


n-1       n- 
2  2 


^     0  /^ 


•  1        n-1 


m^-q     m^-q 


2  2  2 

^2^^     ni^-q 


2 r    ""i/   \  ^ 2"  ■'"2 
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i,      i       I. 


I  l^        I        £. 


(B.6) 


0        0        0/ 


V-^i 


tn,        m       m. 


Now  consider  M'  =  <n-qm    |r  Y |n~qn  >.       This  is  the  same  perturbation  on 

Im 
the  state  which  is  shifted  opposite  to  the  original  state,  by  a  uniform 

electric  "ield. 


m  +m 


m,  — 5—  +q-n      (2X,  +1 )    (21  +1 )    (2I+I ) 

M-    =     I       R,    ,    (^1)    \-\)      ^  /  i-  ^ 

l^l^     ^^2 


4ir 


/n-1         n-1 
/    2  2 


8 


m^+q       m^-q 


n-1        n-1        I 


2\   11^    I    lA  n^     i     i^ 


m^+q     m^-q 


2  2 


Mo    0 


0  /  \  -m         m       m^ 


Recalling  the  symmetry  property 


(B.7) 


m.     m^     m^ 


m„     m,     ra- 


ce.8) 


and  q  is  an  integer  so    (-1)"^*^  =   (-1)"^f  one  obtains 


n;-  : . 
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m. 


m^+m^ 


M'    =     I       R„    „    (-1)    \-1)      ^  "  "/    ^""^ 


-q-n      (2£, +  1)^(2)12+1)^(24+1) 


8.    2. 


4lT 


rn-1  n-1 


2  1 


/- 


n-1       i„\ 


m^-q       m^+q 


2     -^ 


\ 


m^-q     m^+q 


-m. 


2  2        '"2  / 


X  ^  X»  X  ^l       /  **  1  *•  O 


0       0     0/   \ -m         m       m    . 


(-1) 


n-1       n-1        .  n-1       n-*1 

2  2  *1  2  2  2 


(B.9) 


The  last  phase  factor  is      (-1)            .       But/  /      =0  unless 

vo  0     0/ 

\      2  I 

l^   *  l^  +  I     is  even,  so     (-1)              =   (-1)    .  The  result  is 


M'   =    (-D^M. 

If  I  is  odd     («,=1,   dipole;    1=3,   octapole),   states  with  given    |q|   shift 
symmetrically,  either   both  out  or   both  in  towards   the  center. 

If  1  is  even,  states  with  given    |q|   shift  the  same  way   (either  all  to 
the  blue  or  all  to  the  red)  and  so  an  asymmetry  results. 


I 


i  " 


APPENDIX  C 
CALCULATION   OF  Q 


zz 


In  the  spherical  state  representation,  the  matrix  elements  for  Q      , 
within  a  given  principle   quantum  number,   can  be  written 


<n£m|Q^^|nJl'm'>  =  ^2e/ -^  <nim  |r^Y*Q|ni'm'> 


■'^=^rt^-V„,.it/CWr..^»]' 


(C.I) 


Consider  angular  integral,     I    . 


^.  =   ^-^)"'/  VmWrm.^" 


a      «.'      2\     fi      «,'      2 


(-1)™  /   (2£  +  1)(2il'f1)5 

'    -m       m'        0/      I  0       0         0 


X,  £'        2\      /£       i'        2 


(-l)"  /   (2il-^1)(21l'-t-1)5  ^ 

-m       m         0/\0       0         0 

(C.2) 

Note  that   the  result  is   zero  unless  £  +  i'   +  2     is  even  and     (2,,Jl',2) 

satisfy  the  triangle  inequality. 

For  the  m  =  0  case,  the   3-J   symbols  have  a  simple  form. 

!.        H'        2  \  1/0 

0      0       o/'  ^^'^^'- 


(J/2)!  t 

(j/2-S.^)!(j/2-Jl2)!  (J/2-2)!  -.'^'i^ 
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where  J  =  X,  +  i'+2  is  even  and  triangle  inequality  is  satisfied. 

For  the  |m(  =  1  case  and  |m|  =  2  case  (which  is  only  necessary  if 
fine  structure  is  included)  the  3-j  symbols  can  be  calculated  by  the 
method  outlined  in  Vidal,  et  al .  N.B.S.  Monograph  116  (1970).  Now 
consider  the  radial  integral 


I  ^^  ^'Wnl'  • 


0 

We  employ  a  method  'jsed  in  Bethe  and  Sal  peter    (1977)   to  find  the 

normalization  of  the  radial  hydrogenic  wavef unctions  which  are 

proportional   to     e  ^     p  L  ^.    (p)  where  p  =  —  r. 

Using  the  generating  function  for  Laguere  Polynomials   (Abramowltz  and 
Stegum   1972) 

„     ,  ,  (,3)P3-pS/(1-3)  -       lP(p) 

(l-s)  q=p 
one  obtains 


,j.,/2£.k)!(2..k')!    ^^        P  '2£.k^        2'^^k.(P)^P 

r     -p   2U2      (^s)^'^-^e-P^^^^-^^  (,^^2U1^-pt/(1-t) 

=   J      e      p  —2 dp 

0  (l-s)"^^  '^  (1-t)'^'^   '^ 

(st)"  f        2!l+2  -p(1+s/(1-s)+t/(1-t))^ 

°  2F2 2F2  J      P         e  <ip 

(1-s)'^'^  ^(1-t)'^'^   '^     0 

(3t)^'''''^2il->-2)!  (1>-3)(1-t)    ^^'"'^ 

'   (1^3)2i*2(1.t)2^*2  ^^-^^^ 


(st)^^^\..2)i(1-3)(1.t)(^)^^*3 


1-3t' 


fe 
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j=o  •^• 

(st)2^*^1-3)(1-t)      i      (3t)J   lI!:lllilL   .  (c.i*) 


Both  the  left   and  right  hand  side  are  now  power  series   in  s  and  t  and  so 
must   be  equal  term  for  term.      Equate  the  coefficients  of  terms  with 

3  t 


25+k      00                                                                               -                _  /^     n 
(st)^*-  ^    r     -p    2Z+2  ,2l■>■^,    ,2    .           ,    vA2)l+k      (2Z^k->-1)!      ^    ' 
^ — P  J    e     P           L       ^(p)     dp  =    (st)             + 

(2)l+k)!        0  2X,+k  (k-1)!  (     (2il^^k)! 

(k-2)! 
(C.5) 


The  first  terra  on  the  right  hand  side   comes  from     j   =  k-1.       The 

-  -  (28,  +  k)' 

second  term  comes   from     j   =  k-2     and  is  0  if     k  =  1,      ,        . j    ,   otherwise, 

jVpp2^^V^*^p)^dp   =  iii4lll     (2i.k.l)    J  °  .  (C.6) 

Let     n  =  k+i .       Then  [   ]   =  2n  whether  k  =   1   or  k   >   1 . 

J^e    p       L^^^  (p)    dp  =  j^^rpTTT  t^J-  ^^-^^ 


This  determines  the  normalization  factor.     Now  consider 


n 
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22.+l<    29  '+!<'  00 


kk 


"   ^0  iiT2lTior-^  J,,    (2£'.k')!    '  "^P 


22.+2  2{.'+2  P+t'+S 


s^^*^^^'*^^H..4)!  (i-3)^-^^3(,.,).^z-3  I  jf-^;:;:^;;,  ot) 


Case    1)     Z  =  jl' 


Equate  terms  s^t"". 


(C.8) 


mm      * 
3  t      r  ^-p   2t  +  i4. 21+1-    ,,2a  +  1,    ,. 

(m-2Jl-l)!        ^    (m^2Jl-2)!        ^    (in-2Jl-3)!         {m-2l-H)l      '        ^^'^' 


Therefore 


f     ^-p   2a  +  i<,  2]l+l  ,    w2i+1  ,    ,^ 
J      e  ^p  L^        (p)L^        (p)dp 


(^i2l~])l      ^^^"^  '^  ^°^^^  ""    (^8il^-128,  +  i<)m  +    (-81^+1 2£^-J4i)] 


(CIO)  " ''^' 
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Let   m  =  n  +  2,,     put   in  the  normalization,  and  convert  to  an  integral  over 

r  to  get , 

oo  2  a      2 

J     VnJl'^'^'^''  =  |-  [5n^  +   1   -  3l(ilH)]    [^]    .  (C.11) 


Case   2)      £•    =  £   +   2 

The  same  procedure  yields 


a„   2 


/     Rj^j^R^j^^/   dr  =  I  n2/(n^-(il  +  1)^(n^-(il  +  2)^)      (^]      .  (C .  1 2) 


This  is  sufficient   to  give  the  matrix  elements  of  Q_     in  spherical 
representation.     If  it  is  desirable  to  have  Q^^  in  parabolic 
representation,  eqn.   (B.2)   can  be   used. 


APPENDIX  D 
THE  INDEPENDENT  PERTURBER  MICROFIELD 

The  probability  distribution  function  (microf ield)  for  finding  an 

electric  field  of  strength  £  in  a  plasma  is  given  by 

Q(e)  =  <6(£-E)>,  (D.I) 


where  <...>  denotes  a  plasma  average  and     E  ^  i  E.   is  the  sum  of  ion 

i 


fields    (possibly  screened). 


Q(e)   =  lim       1  |...Jd3^   e'^^it-E)  (D.2) 

N,V+« 
N 
V   =  " 


where  Z  is  the  conf igurational  partition  function.  We  will  assume  that 
the  thermodynamic  limit  exists  and  drop  its  explicit  mention  hereafter. 
Inserting  the  integral  definition  of  a  delta  function, 


5(M)  =  -md\ 

(2Tr)^ 


t^t^    _       1         f^3„      ^ii-if-E) 


yields 


Q(^)   =_>        Jd3ie-"-^T(l)  (D.3) 

(2Tr)^ 


where 


■»5f<l 
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Tih   =l/d3^   e-^V^-^ 


(D.H) 


It  is  easy  to  see  that  T(ll)  is  a  function  only  of  the  magnitude 
of  X,  and  not  its  direction.  This  allows  the  angular  integrations  in 
(D.3)   to  be  done  readily,  yielding 


Q(e)    =  -^  j     6i   I  TiDslr.tl. 


(D.5) 


2tt  e      0 


The  Independent  Perturber  approximation  can  now  be  used  to  find  T(il). 
This  approximation  contends   that   interactions  between  the  radiator  and 
perturbing  particles  are  more  important  for  determining  the  microfield 
than  are  interactions   between  different   perturbers.     Therefore  all 
perturber- perturber  interactions  are  neglected,  leaving  only  central 
interations . 


T(l) 


.k 


I  -Sv.+il'E. 
1  1 

3N       i 
re 


h 


I  -ev. 


3N       i 
re 


ilTT  /  dr  r^  e'^''^'^V(£E(r)) 
0  " 


^ir  J  dr 
0 


2     -Bvir) 
r     e 


1   - 


0 


4Tr  /  dr  r^  e'^'^^'^^j    (JlE(r))-1) 


^TT  /  dr 
0 


2     -Bv(r) 
re 


J 


^ 
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where  R   is  the  radius  of  the  system    (assumed  to  be  large).     For  large  r, 
v(r)  ■*■  0  so  the  denominator  is  essentially  the  volume  of  the  system. 

T(o   =  lim  [1   -  ^  /"dr  r^  e'^^^-^^j    (ilE(r))-1)]^ 
N-*-"  0 

00  .         . 

=  expIHirn  J  dr  r^  e'^''^'' ^j    (tE(r  )  )-1 )} . 
0 

This  is     T(J!,)   in  the  Independent  Perturber  Model.     It  is  interesting  to 
note  that     T     (X,)   is  exactly  the  same  as   for  the  first  term  of  Baranger- 
Mozer  formalism  except   that  for  the  I. P.  model   v(r)   is  the  bare  potential 
and  for  B.M.   it  is  screened. 


w- 


APPENDIX  E 
CALCULATION  OF  AN   IMPORTANT   INTEGRAL 


From   (I.B.9),  we  have 


2        2 
sin  k/e   +e.    -2ee.u 


I_  -  JdQ(1^3u   )  ■  .  (E.I) 


2        2 
k/e   +e.    -2tt.\i 


From  Abramowitz  and  Stegun    (1972), 


3xn  k  /e   +e.    -2eG.u  " 

=      I    (2n+1)J    (k£)j    (ke^)P^(u).  (E.2) 

,     /  2       2   „  n=0 


Also  (1-3m^)  =  -2P  (u).   Therefore, 


1 
I  =  -4iT  /  dp  P  (u)  I   (2n+1)j  (ke)j  (ke  )P  (m).  (E.3) 

"       -1        n=0 


Since 

[^   dp  p,.(u)p,(u)  =  2ilT«u-'  ^^''^ 


one  obtains 


I^  =  -8-n3^{kE)3^{kz^).  (E.5) 
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APPENDIX  F 
NUMERICAL  PROCEDURES 

In  this  appendix,  the  numerical   procedures   used  in  calculating 

(I.B.12)  will   be  discussed. 

-ax  2  2 

2^-^l('<)  r         g(x)  ^5j—  (Uax+     ^y-  3^{kz,M) 


<f>     =  ^6e  J     dk  k  e  J2('<e)   J     d 


0  2  0  f'  "^/^^ 

^  ^  J     dk  k  e  sin(ke) 

0 


(F.I) 


Let  W(k)   be  defined  as 


»  -ax  2  2 

W(k)    s  J     dx  g(x) (1    +  ax   +  ^-^  j    (ke(x)).  (F.2) 

0  ^  ^ 


This  integral   is  done  using  trap  rule  integration.     The  outer  integral   is 
broken  into  small   intervals. 


j"  dk  >  I  /'^^^    dk. 
0  i     k. 

On  the  interval    (k^^.k^^^)  let  W(k)   be  approximated  by 

W(k)   -  AW.k     +  BW.k     where  the  constants  AW.,   and  BW.   are  determined  by 
11  J-  1 

2 
matching  the  endpoints  of  the  interval.     (Note  that     W(k)  -  k     for  small   k 

so  the   quadratic  first  term  makes  sense.)     Similarly  let  nh^(k)   =  AH^^k  + 

BHj^  on    (kj^.k^^^).     This  yields 
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-6el  J^       dk  k  e  {[^^  -  — jsin(ke  ) ^-^     }(AW.k+BH.k) 

..^  i        i  k^e^  k  e 

<r>     =  — 


k.^  AH.k+BH. 

I  J  dkke  sin(ke) 

1      '^i 


(F.3) 


Let 


k.  AH.k+BH. 

S.(N)   =  I  ^""^   dk  e     ^         ^   \<~^   3in(k£)     and 

^  k. 

1 

k.  AH.k+BH. 

C.(N)   =   /  ^*'   dk  e     ^         ^   k^" 'cos(ke).  (F.4) 

^  k. 

1 

Then 

AW.  BW.  3AW  3BW.  3AW. 

<f>      =  -6e    I S.(4) i  S.  (3)    +  — T^  S.  (2)    +  — ^  S.  ( 1 ) ^i  C.  (3) 

e  .ei  ei  3       i  3i  2i 

1  e  e  e 


3BW. 


-  C.(2)    /I  S.(2).  (F.5) 


2       i'    '      t-     i 
e  1 


From  Gradshteyn  and  Ryzhik   (1980),  we  have 


AH.x+BH.      N      /    ,sK+1,„    ,,,    N-K  ,k.    , 

_    ,.,■,               1           1      r      (-1)         (N-1)!x               ,     ,             „.  s  i  +  1 

S.  (N)   =  e  2, — j-7Tj~  sin(ex  +  Kt) 

^                                    K=1    (N-K)!(AH.    +e'^)^^  '  k. 


and 


AH.x+BH.      N      ,.,,sK+1  ,„   ,,.    N^K 


fk 


,    ,              i          i     V-      (-'D        (N-1)!x  ,       ^  „.  V 

C.(N)-e  I     _— ___  cos(ex  +  Kt) 

^  K=1    (N-K)!(AH.    +e^)'^^'^  k. 


where 


i  +  1 


(F.6) 
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AH. 
.  .  e  1 
t  =  arcs  in  =  arccos  


/AH.^+c^  /AH.^+e^ 

1  1 

Let 

SL   (N)   =     I     Irrlfr, '  3in(ek     +Kt), 


M  V     1         1/    '^'"^ 

and  similarly  for  CL^(N)  and  CU^(N)  except  use  cos  instead  of  sin.  Then 


AH.k.   +BH.  AH.k.+BH. 

S^(N)  -  (N-1)![e  ^  ^^'   ^SU.(N)  -  e  ^  ^   ^SL^(N)) 


and 


AH.k.   +BH.  AH  k.+BH. 

C.(N)  =  (N-1)!(e  ^  ^■''   ^CU.(N)  -  e  ^  ^   ^CL.(N)).         (F.8) 


SL,   SU ,  CL,  and  CU  are  calculated  as  above;    then  Sj^   and  C^^   are  formed. 
These  are  used  in   (F.5)   to  get  <f>    .       If     ke  <<    1,    it   is  appropriate  to 
use  an  expansion  for     j-Ckc)  and  3in(ke).       In  this  case 


k.    ,  -  AH,k+BH.    ,22  _ 

■-6e  I  j^^*^   dk  k^e     ^         ^(iL£_)(AW.kWk) 

<^>s'  T-| APTk^BH: '- —'  (^-9) 

I     J   ^'^    dk  k  e      ^  ^   ke 

1        k^ 
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Integrals  of  the  following  form  are  needed 


i  +  1  n                    .    r.-r 

f  ^    '  ax+b  n         ax+b     r  ,    ..r  n!x 

I  =  J  dxe         x=e             2.  ^-1)     

X.  r-0               (n-r)!a 

1 


r+1 


1+1 


(F.10) 


These  can  be  easily  coded,   but,   for  small  ax,  the  following  numerical 
problem  still   exists; 


n  n-r 

ax+b    ,      r      /    i\i"   X 

I  =  e         n!      1      (-1)     -p7^ 

r=0  a        (n-r)! 


X.     , 

1  +  1 


ax+b    ,      r      ,   1^n-J  ^ 

=  e         n!      2.      ("1 )         


j=0 


■ .     r.- j  + 1 

j!   a     ^         X. 


X. 


i  +  1 


(-lAl    e^^:^     r      (-ax)^ 


n+1 


J-0 


'i  +  1 


(-1)"n!    e^'^:^     ^-ax 


n+1 


-     I 


(-ax)J 


1  ! 


j=n+1 


i  +  1 


(-1)  n!    e 
n+1 


1   ^  e 


ax 


I 
j'^n+l 


(-^ax)- 
J! 


X. 


i  +  1 


(F.11) 


For  small  ax,   the   1   dominates  the  summation;    but   1    is  just  a  constant 
and  so  drops  out  when  taking  the  difference  between  upper  and  lower 
limits,  leaving  very  little  accuracy.     One  may  approximate 


T        ("D       n!      ax+b     y  (-ax)^    I 


n+ 1  '•  i '         '  X 


.      n+1           n+2 
ax+bj-jc ax i 

^  n+1   '  (n+1)  (n+2)  *   '"i 


i  +  1 


X. 

1 


n+1                             2  2  /'^i  +  l 

< r.       ax         a   X 1 

n+1    ^'   "  n+2  ■*■   (n+2)(n+3)   ^   ••'M  *  (F.12) 


x. 

1 

This  is  a  better  expression  to  use  for  computation  when  ax  <<   1 . 

For  large  values  of  e,   the  calculation  of  <f>     detailed  above  has 

e 

non-physical  oscillations.     This  is  a  common  feature  of  numerical  Fourier 
transforms.     To  avoid  this   problem,  we   use  an  asymptotic  form.     For 
large  e     values,  the  nearest  neighbor  limit   (<f>         =  -2e^^^/>/z)   is 

valid.     Divide  the  expression  for   <f>     by  the  nearest  neighbor  form  to  get 

1  /2 
a  function  that  goes   to   1   for  large  e   and  goes  as  e         for  small  e.       This 

can  be  fit   by  a   Fade  approximate  of  the  form 


'''^  A/F  +   BE   +   C-3/2 

'NN 


<f> ^  ;= 372  ^^-^3) 


where  A,   B,   C,  D,  and  E  are  fitting  parameters.     These  parameters  are 
determined  by  minimizing  the  function 


^-^t  A/e   +  Be  +  Ce^^^ 


''>^  1    +  D/F  +   Ee  +   Ce3/2 


^NN 


for  small  values  of  e .       A  code  "Simplex"   provided  by  Dr.   Robert  Col dwell 
does  the  minimization. 


'-:/^ 


APPENDIX  G 
ELECTRIC  FIELD  BEHAVIOR   OF  THE  ENERGY  LEVELS 
OF   HELIUM-LIKE   IONS 

From  equation   (IV.A.1)  we  have 


J((D,e)    >  -  ;^  Im  Tr  DJoj-CH^-oiq)   -  p-£   -  VCw)}"^    .  (G.I) 

a 


Using   (IV. A. 2),  we  can  rewrite  this,   for  a     Ly-a  line,  as 


J(a),e)    = Im  Tr  dIaoj  +  u     6(!!.,0)   -  P'e   -   ip} 


^         a  ^P 


where 


Ao)   -  a)-(Z-1)^(1-  i^)   -  a)(2,1,Z) 
n 


u)        =  (ij(2,1,Z)   -  a)(2,0,Z) 
sp 


(G.2) 


ir   -  ?/(a)). 

We  have  expressed  ^(to)     as   ir     to  show  that  it  is  pure  imaginary  and  that 
since  it  is  essentially  constant  out  to  the  electron  plasma  frequency,  we 
will  treat  it  as  a  constant  here.     We  have  also  chosen  the  isolated  p 
state     (i-l)   as  our  energy  reference  point.     We  will  use  the  spherical 
representation,    |n£ra>,   and  chose  the  states,   for  our  matrix,   in  the 
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fol lowing  order    |200>,    |210>,    |2n>,   and   |21-1>.      In  this  representation, 

the  resolvent  R,    (aw  +  o)     6(Jl,0)  -  p«e  -  if}      ,     has  the  following  form 

sp 

(Smith  1966) 


Au-o) - 

sp 

-if 

E 

0 

0 

E 

A(D-ir 

0 

0 

0 

0 

Aio-ir 

0 

0 

0 

0 

Auj-ir 

-1 


(G.3) 


Mu-ir 


-E 


where 


Aw+o)     -ir 
-E  ^^—  0 


X  -   (Au+(D     -ir)(Au)-ir)-E     and 
sp 


E  -  3eaQe/h, 


AwQ-ir 


0 
1 


Au-ir 


K  >  • 


The  matrix  D  has  the  form 
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0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 
0 


(G.4) 


This  yields 


J((u,e) 


=  -  1  lm{ 


AiD+o)     -ir 
sp 


A(D 


%l 


(G.5) 


The  second  term  corresponds  to  the  unshifted  components.     Since  we 
are  interested  in  the  field  dependence  of  the  shifted  components,  consider 
only  the  first  term.     The  peaks  in  JCu.e)   corresponding  to  the  energy 
levels  of  the  shifted  components  occur  at   points  for  which  the  real   part 
of  X  is  equal  to  zero.     Neglecting     r,     these  are  given  by 


2  2 

Adi   +Au  0)       -  E     »  0 
sp 


Ao) 


2        2 

sp       sp 


(G.6) 


In  order  to  better  understand  this  expression,  we  shall  consider  two 

limits.  First,  for  u   >>  E, 

sp 


ft. 


*»■.,  ■' 


-111- 


Au  '(  .  (G.7) 

i_ 


For  0)   <<  E,   we  have 
sp 

2 

Ao)  «<  .  (G.8) 

01       0) 

__  _  sp  _  sp 
*"    2     8E 
The  first  limit  is  the  small  ion  field  limit.  Here  we  see  the  states 

2 
shift  as  z    .       The  second  limit  is  for  large  ion  fields  and  shows  linear 

dependence.  The  general  expression,  (G.6),  goes  smoothly  from  the 

quadratic  Stark  effect  for  small  fields  to  the  linear  Stark  effect  for 

large  fields. 
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